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A natural extension to the standard model of particle 
physics is a fourth generation of quarks and leptons. The 
inclusion of a fourth generation provides a source of CP 
violation in decays and can accommodate a heavy 
Higgs boson fT| [51 . Searches for fourth generation quarks 
at the Fermilab Tevatron have constrained the mass of 
up-type quarks (W4), that decay as U4 — > Wq, where q is 
a generic down- type quark, to be m„_j > 340 GeV/c^ at 
95% confidence level (CL) [3], while limits on the mass of 
down-type quarks (d^ decaying via d^^ — > Wt are 771^4 > 
372 GeV/c^ at 95% CL 4J. 

Following the trend of mass hierarchy in the standard 
model, the least massive and therefore most accessible 
particle of this fourth generation may be the neutrino. 
Such a neutrino need not be solely a Dirac or Majorana 
state, but may be a mixture of the two [51. This leads to 
two mass eigenstates A^i and A'2 , where A'2 is the unstable 
heavy eigenstate and A^i is the stable and least massive 
eigenstate of the fourth generation neutrinos. These par- 
ticles would partially evade the neutrino mass constraints 
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from Z width studies at LEP [3]. 

The dominant production mechanism of iVi would 
be via a Dreh-Yan process, -pp — ■^/7* — > N2N2 — > 
NiZNiZ, giving a final state of two Z bosons and large 
missing transverse momentum. This signature is shared 
by several other interesting new physics processes, most 
notably supersymmetric production, X2X2 ~^ ^Xi^Xi^ 
where Xi ^^'^ X2 ^-^^ neutralinos. We consider the 
mode in which one Z decays hadronically and the other 
decays leptonically, giving a detector signature of two 
charged leptons, two jets and large missing transverse 
momentum. For this search we use pp collisions at 
= 1.96 TeV corresponding to 4.2 fb~^ of integrated 
luminosity collected by the CDF II detector. 

Events were recorded by CDF II [3 [S], a general- 
purpose detector designed to study collisions at the Fer- 
milab Tevatron pp collider. The CDF II detector is com- 
posed of a charged-particle tracking system immersed 
in a 1.4 T magnetic field consisting of a silicon mi- 
crostrip tracker and a drift chamber. Electromagnetic 
and hadronic calorimeters surround the tracking system 
and measure particle energies. Drift chambers located 
outside the calorimeters detect muons. 

The data acquisition system is triggered by e or can- 
didates with transverse momentum px, greater than 18 
GeV/c. We retain electron and muon candidates with 
pseudorapidity [5] I77I < 1.1, pt > 20 GeV/c and that 
satisfy the standard CDF identification requirements [S] . 
For muons, the track fit per degree of freedom is used 
to reject poorly fit tracks likely resulting from charged 
pion and kaon decays in flight. Electrons from photon 
conversions are suppressed by rejecting electron candi- 
dates with a nearly coUinear intersecting reconstructed 
track. Jets are reconstructed in the calorimeter using 
the JETCLU [TD] algorithm with a clustering radius of 
0.4 in azimuth-pseudorapidity space. Measured jet ener- 
gies are corrected to account for 77-dependent variations 
in detector response, calorimeter coverage, and the ex- 
pected contribution from additional pp interactions in the 
same event [TT]. Jets are selected if they have pr > 15 
GeV/c and |r?| < 2.4. Missing transverse energy p^, ^t, 
is reconstructed using calorimeter and muon information 
including the corrections described above. 

To isolate the ZZ signature, we require two opposite- 
charge, same-flavor lepton candidates (e or fi) with pr > 
20 GeV/c for which the lepton-pair invariant mass is 
consistent with decay from a Z boson: to« £ [76, 106] 
GeV/c^. Additionally, we require at least two jets, each 
with pt > 15 GeV/c and {rjl < 2.4, and without identified 
secondary vertices resulting from 6-hadron decay [13] . 
The ZZ + signature has the further requirement of 
large ^t, varying with hypothetical A'^i and N2 masses, 
as shown in Table HTl 

The dominant background in the resulting sample is 
production of a Z boson in association with two jets from 
initial state radiation. We model this background using 
ALPGEN [13] to describe the hard process and pythia [15] 
for the showering and hadronization. This background is 



strongly suppressed in events with large missing trans- 
verse momentum, as shown in Figure [l]and Tabled and 
is distinguished from the signal by the lack of a resonance 
in the dijct mass, rrijj. 

The second largest expected background is due to W 
boson production in association with three jets from 
initial state radiation, where one jet is wrongly recon- 
structed as a lepton. We model this using an independent 
sample of events containing jets likely to mimic leptons, 
following Ref. fTBl. Additional backgrounds result from 
standard-model production of two gauge bosons, includ- 
ing ZZ, WW, and WZ, as well as tt WbWb, which 
are all modeled using pythia. 
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Missing Transverse Momentum [GeV/cJ 

FIG. 1: Distribution of missing transverse momentum in 
events with the ZZ signature, for expected backgrounds and 
observed data. 



TABLE I: Expected number of events for each source of back- 
ground to the ZZ -)■ e+rjj and ZZ + X1X1 l+Tjj+pT 
signatures, as well as the observed event yield in data with 
4.2 fb~^ of integrated luminosity. The threshold in is 
optimized as a function of the Ni,N2 masses; one example 
(A^i = 125 GeV/c^^'2 = 225 GeV/c^) is shown here. Uncer- 
tainties shown include both systematic and statistical uncer- 
tainty added in quadrature. 



jj and 



Process 


£+£-jj 


> 36 GeV 


WW 


4.4 ± 1.3 


2.7 ± 0.8 


tt 


14.8 ± 3.0 


11.6 ± 2.3 


H^-l-jets 


36.1 ± 16.7 


21.7 ± 12.6 


ZZ 


99.4 ± 20.5 


4.2 ± 0.9 


WZ 


105.6 ± 22.1 


5.2 ± 1.1 


Z+iets 


10171 ± 4422 


94.6 ± 38.5 


Total 


10432 ± 4485 


140.0 ± 40.6 


Data 


10199 


152 



To isolate the double- resonance nature of the ZZ + 
signature, we calculate the distance from the Z boson re- 
constructed mass in the mu—rrijj mass plane, accounting 
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for the relative difference in the resolutions between the 
leptons and jets as well as the observed bias in recon- 



structed 



using the variable 



' mu - mz- 



9n 



(1) 



where mu{nijj) is the reconstructed lepton (jet) pair 
mass, compared to the reference mz-^u — 91.6 GeV/c^ 
{mz^jj = 85.3 GeV/c^) found in simulated events. To 
account for the superior lepton resolution, the dilepton 
and dijet mass differences are scaled by factors related 
to the resolutions: gu = 10 GeV/c^, gjj = 15 GeV/c^. 
The uncertainties of these reference values are small, and 
may be neglected. The distribution of Am for data and 
simulated background and signal is shown in Figure [2] 

We model the production of the N2 signal and its 
subsequent decay into A^i over a grid of masses in 
the {Mni, Mpf2) plane using madgraph |T7j with the 
CTEQSl [18] parton distribution fuctions; pythia [15] 
is used for the showering and hadronization. To sup- 
press the large backgrounds expected from standard- 
model sources we require large ^t', as the expected mag- 
nitude of missing transverse momentum depends strongly 
on M^ri and M7V2, we vary the selection threshold of 
to optimize for sensitivity at each {M^i, M]\f2) pair con- 
sidered, as seen in Table [TT] The acceptance for each mass 
point can be seen in Figure |3] For each point in the mass 
grid, we form template histograms as a function of Am 
for the expected signal and background, as displayed in 
Figure [2] 



M^i = 225 GeV/c" 
M^^ = 125 GeV/c- 

^1 Background 

I I Signal X 10' 

Background Uncertainty 

• Data 



A m 



FIG. 2: Distribution of the variable Am, defined in the text, 
for expected background, observed data and an example sig- 
nal (scaled by 10*) in data with 4.2 fb~^ of integrated lumi- 
nosity. This example uses a missing transverse momentum 
threshold of > 36 GeV, optimized for this {Mni, Mn2) 
mass point; see Table [ll] Background uncertainties are statis- 
tical and systematic added in quadrature. 



In addition to the templates formed for the nominal 
expectation, we form alternate templates that incorpo- 
rate the effects of systematic uncertainties under ±lcr 




220 240 260 2S0 300 320 



M„ [GeV/c-] 



FIG. 3: Acceptance of the ZZ -f signature, including 
3K{ZZ — >■ Uqq), as a function of the masses of the fourth 
generation neutrinos, A'^i and A^2. The threshold in is 
optimized at each point on a grid in this plane. Linear inter- 
polation is performed between the grid points. The apparent 
structure in the plot results from statistical fluctuation. 



variation. Fitting to these templates using the maxi- 
mum likelihood method, we extract the best-fit signal 
cross section, (TAr2. Systematic uncertainties affecting the 
shapes of templates, including uncertainty in the jet en- 
ergy scale [TT], QCD radiation, PDFs, (square of mo- 
mentum transfer in the interaction) and uncertainty in 
lepton energy resolution, are accounted for as nuisance 
parameters in our likelihood. The dominant source of 
systematic uncertainty in this analysis is uncertainty in 
the jet energy scale (40%), which can significantly modify 
the number of jets in background processes that pass the 
Pt threshold, the location of the rrijj resonance in the sig- 
nal process, and the measured in an event. The sec- 
ond largest systematic uncertainty is due to uncertainty 
on the theoretical normalization of the background rates. 
Finally, we apply the unified ordering principle jl9| for 
the Neyman construction to create confidence intervals 
in the true value of 0^2 for each 7V2, A'^i mass point. 

We find the candidate events in the data to be con- 
sistent with expected standard-model backgrounds and 
thus set upper limits at 95% GL on the cross section for 
pp N2N2 N\ZN\Z. Theoretical cross sections for 
each mass point are presented in Table |ll] along with 
their respective expected and observed limits in our data 
sample. The expected and observed cross section limits 
can be seen in Figure |4| and Table [TT] 

In summary, we have performed the first search for new 
phenomena in events with two reconstructed Z bosons 
and large missing transverse momentum. This signature 
is sensitive to processes pp — ^ X2X2 — > ZZX^X^-, where 
X2 is an unstable particle decaying as X2 — > ZX\ and 
X\ being undetected. The particles X\ and X2 may be, 
among other possibilities, fourth generation neutrinos or 
supersymmetric particles. A specific model in which X2 
and X\ are fourth-generation neutrinos is used without 
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TABLE II: Acceptance of the ZZ + ]pT selection for vary- 
ing thresholds in ]pT optimized for each point in the Mni, 
Mn2 mass plane. Also shown are the median expected and 
observed 95% CL upper limits on the cross section (o-jv2) in 
data with 4.2 fb~^ of integrated luminosity, as well as the 
theoretical prediction |17l I20| . 



Mjvi, Mjv2 Cut Acceptance (Tjv2 [fb] 



[GeV/c^] 


[GoV] 


[%] 


Theory 


Exp. /Obs. Limit 


75, 175 


37 


0.99 


0.51 


511 / 702 


75, 200 


68 


1.02 


0.21 


292 / 369 


125, 225 


36 


0.85 


0.16 


684 / 1088 


75, 225 


92 


0.93 


0.081 


156 / 273 


75, 275 


118 


1.01 


0.015 


94 / 132 


125, 300 


119 


1.06 


0.013 


99 / 138 


175, 300 


80 


0.96 


0.022 


171 / 315 


125, 350 


156 


1.05 


0.003 


75 / 48 


225, 350 


80 


1.05 


0.006 


190 / 297 


75, 350 


167 


1.06 


0.001 


71 / 55 



loss of generality. In the final state in which one Z bo- 
son decays to two charged leptons and the second de- 
cays hadronically, we find agreement between the data 
and the standard-model expectation using data from 
proton-antiproton collisions with 4.2 fb~^ of integrated 
luminosity. Based on the results in Table we report 
95% CL upper limits on the cross section of the process 
pp X2X2 — > ZZXiXi ranging from 50 fb to 1 pb 
depending on the masses of Xi and X2 ■ 




M^j [GeV/c-] 



FIG. 4: Upper limit at 95% CL on the cross section of 
pp iV2A^2 NiZN-iZ in data with 4.2 fb"^ of integrated 
luminosity as a function of the masses of A''i and N2. Top 
shows median expected limits; bottom shows observed limits; 
see Table HH 
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